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We use a neutron spin-echo method with yueV resolution to determine the life- 
times of spin waves in the prototypical antiferromagnet MnF2 over the entire 
Brillouin zone. A theory based on the interaction of magnons with longitudi- 
nal spin fluctuations provides an excellent, parameter-free description of the 
data, except at the lowest momenta and temperatures. This is surprising, given 
the prominence of alternative theories based on magnon-magnon interactions 
in the literature. The results and technique open up a new avenue for the in- 
vestigation of fundamental concepts in magnetism. The technique also allows 
measurement of the lifetimes of other elementary excitations (such as lattice 
vibrations) throughout the Brillouin zone. 
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The concept of elementary excitations is one of the basic pillars of the theory of solids. 
In the low-temperature, long-wavelength limit, such excitations do not interact and have an 
infinite lifetime. For nonzero temperatures and momenta, the lifetimes of elementary excita- 
tions are generally limited by collisions with other excitations, with important consequences 
for the macroscopic properties of solids. For instance, the thermal expansion of solids can 
be understood as a consequence of collisions between lattice vibrations (phonons). Because 
of their comparatively simple Hamiltonians, magnetically-ordered states are excellent test- 
ing grounds for theories of elementary excitations and their interactions. Despite this, the 
damping of spin waves (magnons) in antiferromagnets has remained an open problem for 
four decades. Theoretical calculations of magnon lifetimes have been carried out since the 
1960's, with intensive development occurring on several fronts in the early 1970's. However, 
these activities ground to a halt by the mid-1970's due to the lack of appropriate experimental 
data, namely, from momentum-resolved measurements with sufficient energy resolution. The 
only low-temperature data available were taken with q = 0, in antiferromagnetic resonance 
(AFMR) and parallel pumping measurements ID El- Because of the limited energy resolution, 
momentum-resolved data from neutron spectroscopy |l3l, on the other hand, were confined to 
the critical regime extremely close to the Neel temperature (Tn), where most theories do not 
apply. Until recently, no other experimental techniques were available which permitted high- 
resolution measurements of excitation lifetimes at low temperatures over the whole Brillouin 
zone. We report on a new neutron spectroscopy method with yueV resolution which is used to 
measure spin wave (magnon) lifetimes in the prototypical antiferromagnet MnF2 over the tem- 
perature range from 0.04 - 0.6 Tn. The results subject long-standing theoretical predictions to 
a first experimental test and hold promise as a novel probe of elementary excitations in quan- 
tum magnets. The technique is also widely applicable to other elementary excitations such as 
phonons and crystal field excitations. 
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The determination of magnon lifetimes at low temperatures requires an energy resolution in 
the jjeV range, about two orders of magnitude better than that achievable by standard neutron 
triple-axis spectroscopy (TAS). We have obtained the requisite gain in resolution by manip- 
ulating the Larmor phase of the neutron spin with magnetic fields. The TRISP spectrometer 
(Fig. lA) weds the capability of TAS of accessing collective excitations throughout the Bril- 
louin zone to the extremely high energy resolution of neutron spin-echo spectroscopy |l5|. As 
in typical spin-polarized triple-axis spectrometry, the neutrons impinging upon the sample are 
polarized, and the polarization of neutrons scattered from the sample is measured. On TRISP, 
this is accomplished through the use of a polarizing neutron guide and a transmission polarizer, 
respectively. However, in analogy to neutron spin-echo spectrometry, the TRISP spectrome- 
ter also includes regions of effectively constant magnetic field which are produced by pairs of 
radio-frequency (RF) resonance coils inserted symmetrically 1) between the monochromator 
and sample and 2) between the sample and the analyzer [6|. The RF frequencies in the coils are 
tuned such that each detected neutron that creates an excitation lying on the magnon dispersion 
curve has the same net Larmor phase after traversing the two spin-echo arms, independent of 
small variations in the wave vector of the excitation. The neutron spin polarization determined 
at the detector is then a measure of the linewidth (inverse lifetime) of the magnon. In this way, 
the measured linewidth is decoupled (to first order) from the spread in energy of the neutrons 
incident on the sample, which is responsible for the instrumental resolution in TAS. (For a 
detailed description of the technique, see the Materials and Methods section.) 

We chose the antiferromagnet MnF2 for the experiment, as its magnetic ground state and 
excitations have been investigated extensively. The lattice structure and magnetic ordering of 
MnF2 are shown in Fig. IB. MnF2 has the body-centered tetragonal structure, with a = b = 
4.8736 A and c = 3.2998 A. The Mn^+ ions have spin S = 5/2, and the spin in the center of the 
unit cell is oriented antiparallel to those at the comers. The strongest magnetic interaction is 
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between second-nearest-neighbor Mn^+ spins (comer and center spins) and is antiferromagnetic 
0. A weaker, ferromagnetic interaction exists between nearest-neighbor spins (along the c- 
axis). A relatively strong uniaxial anisotropy which is predominantly the result of dipole-dipole 
interactions |l8l|9| causes the spins to align along the c-axis. Tn is 67.6 K. The slope of the 
magnon dispersion is required to set the tilt angles of the RF coils and to determine the non- 
intrinsic contribution to the data (see Materials and Methods). During the experiment, the spin- 
wave dispersion was therefore measured at each temperature at which linewidth data was taken; 
a partial data set is shown in Fig. IC. 

Fig. 2 shows raw polarization data as a function of the spin-echo time r. The spin-echo 
time is proportional to the frequency in the RF coils and the distance between the coils, and 
also depends on the neutron wavelength. In a neutron spin-echo experiment, the dependence of 
the measured polarization on r corresponds to the Fourier transform of the scattering function 
as a function of energy. The data in Fig. 2 are described well by an exponential decay, which 
indicates that the spectral function which characterizes the magnon linewidth is a Lorentzian in 
energy. The difference in linewidth (half- width at half maximum, or HWHM) between the up- 
per two and lower two data sets is in each case only ~3 /ieV, but it can be resolved clearly. The 
upper pair of data sets represents a difference in q of 0.05 r.l.u. at 15 K ifTOll . For comparison, 
the HWHM of the corresponding TAS scans of the lower two magnons in Fig. 2, taken with 
fixed final neutron wave vector kf=l.l A^^, is approximately 100 yueV. 

The raw data were then corrected for instrumental and non-intrinsic effects ifTTl . Figs. 3 
and 4 show the intrinsic magnon linewidth as a function of momentum q and temperature T, 
respectively. The linewidth generally increases with increasing q and T, due to the increasing 
likelihood of collisions with other excitations. However, Fig. 3 also shows that the linewidth 
deviates from this general trend and exhibits peaks as a function of q close to the center and 
the boundary of the antiferromagnetic Brillouin zone. The low-g peak is already present at 3 
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K, the lowest temperature covered by this experiment, and it evolves weakly with increasing 
temperature. This behavior is not described by the dominant magnon relaxation mechanisms 
for which quantitative predictions are available; possible origins will be discussed below. In 
order to facilitate comparison with these predictions, we have treated the 3 K linewidth data as 
a temperature-independent contribution and subtracted it from the higher-temperature data. The 
results are shown in the main panels of Figs. 3 and 4. 

The intrinsic relaxation channel for magnons that has received by far the most attention in 
the literature is magnon-magnon scattering. In an "n-magnon" scattering event, a magnon (here, 
one excited by an incoming neutron) scatters off {n/2 — 1) thermally excited magnons, produc- 
ing n/2 scattered magnons which are in thermal equilibrium with the sample. In the absence of 
defects and external magnetic fields, the lowest-order interaction which limits the magnon life- 
time is 4-magnon scattering. Unfortunately, a comprehensive survey of the literature revealed 
very few theoretical predictions appropriate for comparison with our data, despite the existence 
of considerable work on 4- and 6-magnon interactions. This is either because the calculations 
employed approximations valid only in high magnetic fields (for the purpose of comparison 
with AFMR data), or because strict inequalities that define the range of applicability of the the- 
oretical results are extremely difficult to satisfy experimentally. An analytical expression was 
given by Harris et al. fV2\, who evaluated the contribution to the linewidth from 4-magnon 
scattering processes analytically for the case of single-ion anisotropy with g = lfT3ll . The cor- 
responding result is shown in the bottom trace of Fig. 4. At low temperatures, the temperature 
dependence of the data is considerably weaker than that predicted by this theoretical result for 
g = 0. The best agreement of the magnitude occurs at 40 K, where the experimental result is 
~ 30% larger than the theoretical. Predictions for an anisotropy gap of dipolar origin, which 
would be more appropriate for MnF2, are not available. 

An additional relaxation channel, in which magnons are scattered by thermally-excited lon- 
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gitudinal spin fluctuations, was considered by Stinchcombe and coworkers lfT4lfT5ll . The curves 
in Figs. 3 and 4 are based on this mechanism. (For g = 0, where the contribution of this re- 
laxation is identically zero, we have shown the prediction of the 4-magnon relaxation model, as 
discussed above.) For the larger-g data, the linewidth far from Tn is given approximately by 

rJHWHM) = f^^'',^' , (1) 

where is the magnon energy, q* = 27iq/a, /i* = 2.969 A^, p* = 5.864 A^, and P = l/keT, with 
ke the Boltzmann constant fTS]. The anisotropy parameter a is equal to 0.0184. The exchange 
parameter J(0) = 6.02 meV includes both first- and second-nearest neighbor exchange interac- 
tions. The parameters Rq and Rq, which are both temperature-dependent, can each be evaluated 
using either experimental data or results from mean-field theory ifTSl . leading to considerable 
differences in the magnitude of the calculated linewidth and in its variation with temperature. In 
determining Rq, we used experimental data for the staggered magnetization lfT6l . Calculation 
of Rq from the Brillouin function produces linewidth values which agree at the lowest tempera- 
tures and begin to deviate with increasing temperature: at 40 K, the calculated linewidth is 1 1% 
smaller. For Rq, we used the derivative of the Brillouin function. Calculation of Rq instead from 
experimental data for the parallel magnetic susceptibility ifTTI produces linewidth results which 
are 40% larger at 15 K and 30% smaller at 40 K. 

Given the prominence of the magnon-magnon scattering channel in the literature, the ex- 
cellent agreement between this model calculation and the experimental data is surprising. The 
dominance of the relaxation by longitudinal fluctuations is, however, consistent with arguments 
by Reinecke and Stinchcombe, who estimated that the contribution to the linewidth from 4- 
magnon scattering is only l/z of the magnitude of the above term ifTSll . Here, z is the number 
of neighbors which experience the strongest exchange interaction; 2; = 8 for MnF2, for which 
case z is the number of next-nearest neighbors lfT9ll . As the analytical expression on which the 
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curves in Figs. 3 and 4 are based is valid only at low q, deviation from the data at larger q is 
not unexpected. The general expression for the linewidth resulting from scattering by longitu- 
dinal spin fluctuations lfT4lfT5l should be evaluated numerically at high q to see if the peak as a 
function of q can be reproduced ll33ll34l . 

An explanation of the peak centered at g ~ 0.1 r.l.u. (inset in Fig. 3) requires a differ- 
ent mechanism [20|. An additional potential source of linewidth is the hyperfine interaction, 
which gives rise to the scattering of electronic magnons from nuclear spin fluctuations ll2Tl . 
The contribution from the hyperfine interaction would only be weakly temperature dependent, 
because the nuclear spin system is already highly disordered thermally at 3 K. 4-magnon scat- 
tering terms in which one electronic and one nuclear magnon interact have indeed been shown 
to generate maxima in the linewidth at nonzero q, but estimates of the amplitude of this contri- 
bution are at least an order of magnitude smaller than the observed effect ll2Tl . The crossing of 
magnon and TA phonon modes at g = 0.04 r.l.u. may also contribute to the peak ll22l l23l l24l . 
An additional relaxation mechanism which must be considered as a possible source of linewidth 
is that of magnon-phonon scattering. Experimental estimates of the linewidth due to magnon- 
phonon relaxation in MnF2 in zero field range over three orders of magnitude, but again appear 
too small to explain the observed peak Il25ll26ll27l . A theoretical estimate of the spin-lattice 
relaxation time (which should be of the same order of magnitude as the magnon-phonon relax- 
ation times) corresponds to a linewidth of ~0.5 fieV at 25 K in MnF2 ll28ll29l . Li this theory, 
the magnon-phonon interaction arises from the phonon modulation of the exchange interaction, 
and is dominated by 2-magnon-l -phonon processes. The result varies as T^, which corresponds 
to a linewidth of 1 fieV at 30 K and 5 fieV at 40 K. The maximum potential contribution to our 
data would then be ~ 60% of the linewidth at g = and 40 K. Other mechanisms which may 
contribute to the presence of this peak include 2- and 3-magnon non-momentum-conserving 
processes which originate from scattering from defects 1301131112111 . The linewidth originating 
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from the latter process is peaked at intermediate q. Using parameters derived from comparison 
with data on RbMnFs, its contribution in MnF2 can be estimated to be two orders of magnitude 
smaller than the data OTl . 

The challenge to theory posed by the temperature- and momentum-dependent peaks in the 
magnon linewidth in MnF2 should stimulate new activity in the field of spin wave decay mecha- 
nisms. High-resolution lifetime measurements over the full Brillouin zone in a relatively simple 
antiferromagnet such as MnF2 permit detailed evaluation of proposed processes, which should 
provide a basis for addressing such interactions in more complex magnetic systems. 
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Figure 1 : (A) A diagram of the spectrometer TRISP at the FRM-II. G denotes the polarizing 
guide and AP the transmission polarizer; M and A are the monochromator and analyzer, as in 
TAS. S is the sample and D the detector; VS indicates the velocity selector. The resonance 
coil pairs (Ci and C2) are shown in red, and the mu-metal shielding boxes which enclose them 
in gray. The blue ray represents the path of the neutrons through the spectrometer, from left 
to right on the diagram. (B) The crystal and magnetic structure of MnF2. The gray (smaller) 
spheres represent Mn^+ ions and the green (larger) spheres the F^ ions. The arrows indicate the 
relative directions of the Mn^+ spins on the respective sublattices. (C) The magnon dispersion 
along the qc direction at three selected temperatures at and below 40 K. The data was taken 
on TRISP during the course of the linewidth measurements. The curves show the results of 
fits based on the same spin-wave result used by Okazaki et al. [7|, in which the anisotropy is 
expressed by a single-ion form, and in which the interactions of up to third-nearest neighbors 
are taken into account. 
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Figure 2: Raw polarization data taken at (uppermost data) q = 0.2 r.l.u. and T = 15 K, (middle) 
g = 0.15 r.l.u. and T = 15 K, and (lowest) q = 0.2 r.l.u. and T = 20 K. The lines are exponential 
fits to the data. The corresponding Lorentzian magnon linewidths (HWHM) are 12.4 ±0.8 fieV, 
15.6 ± 0.9 fieV, and 18.6 ± 1.0 /xeV, respectively. 
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Figure 3: Intrinsic magnon linewidth LW at temperatures ranging from 15 - 40 K, as a function 
of q. We have plotted (LW(T, g)-LW(3 K, g)), where LW(3 K, q) is given in the inset (see text). 
The curves show theoretical expressions from Refs. lfT?l and llT51l (see text). Two different 
theoretical expressions are valid for the small-g case, depending on the magnitude of q relative 
to the anisotropy energy. Both expressions apply only to small g; Stinchcombe and Reinecke 
have applied one of them to data extending up to g = 0.2 r.l.u. for MnF2 near Tn 1 15J . However, 
except at g = 0, the theory provides an excellent fit to the data as a function of g up to g ~ 0.35 
r.l.u., both in the magnitude and in the g-dependence. 
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Figure 4: Intrinsic magnon linewidth with < g < 0.35 r.l.u. and q = 0.5 r.l.u., shown as a 
function of temperature. As in Fig. 3, data taken at 3 K has been subtracted. For 0.1 < g < 0.35 
r.l.u., the curves show theoretical results from Ref. [15|. For g = 0, results from Harris et al. 
are plotted [.12J . 
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